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ness modes,  and  the  modal  temperature  coefficients. 
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INTRODUCTION 


Precision  frequency  control  requirements  for  digital  communication  and 
position  location  systems  currently  undergoing  development  make  it  imperative 
that  crystal  resonator  performance  be  improved  in  a number  of  aspects. 
Accordingly,  the  potential  of  doubly  rotated  quartz  cuts  has  begun  to  be 
explored. 3*  For  cuts  on  the  upper  zero  temperature  coefficient  locus  in 
general  (0  * +34°),  and  in  the  neighborhood  of  the  SC-cut  in  particular 
(♦  « 21.9°,  0 “ +33.9°),  a variety  of  effects  having  their  bases  in  nonlinear 
elasticity  have  been  shown,  or  are  predicted,  to  be  reduced  below  the  corres- 
ponding AT-cut  values.  In  addition,  the  static  frequency- temperature 
behavior  shows  some  improvement. 


Among  the  nonlinear  effects  of  interest  are: 

• Force-frequency*  ♦ **»  13-25* 

• Acceleration-  frequency^*  2 7* 

• Resonance  amplitude-f requency*0 » 28-31* 

• Intermodulation*2-35* 

• Mode  coupling-activity  dips-^-43* 

• Dynamic  thermal-f requency^ ,4 , 7-9 

• Film  stress-frequency* *5. 44* 

Additional  references  and  discussion  may  be  found  in  References  6 and  45.* 


Some  of  these  nonlinear  effects  have  received,  or  are  receiving,  theoretical 
and/or  experimental  treatments;  this  technical  report  is  principally  con- 
cerned with  the  force-frequency  effect  which  has  thus  far  not  been  investiga- 
ted in  any  detail  for  doubly  rotated  quartz  plate  vibrators.  This  effect 
relates  the  initial  stress  produced  by  the  mounting  supports  to  resonance 
frequency  changes;  it  contributes  to  long-term  aging  and  is  also  related  to 
the  frequency  excursions  produced  in  shock  and  vibration  environments. 


In-plane  diametric  forces  applied  to  the  periphery  of  vibrating  plates 
produce  frequency  changes  (order  10~?  per  gram)  that  depend  upon  the  azimuth 
angle  in  the  plane  of  the  plate.  If  ip  is  measured  from  the  axis,  then 
it  is  found  experimentally** 11*14-18,20  that  for  the  AT-cut  the  effect  is  zero 
at  values  of  approximately  60°  and  120®.  For  the  IT-cut  at  $ ■ 19.1°, 
Ballato*  found  the  zeros  to  occur  at  i}>  ■ 85°  and  163®  with  a maximum  value 
about  one-third  that  of  the  AT-cut.  This  points  to  a reduced  coefficient  at 
the  SC-cut  as  well. 


In  this  report  we  extend  the  force-frequency  effect  measurement  to  doubly 
rotated  quartz  plates  on  the  upper  zero  temperature  coefficient  locus,  con- 
centrating on  the  SC-  and  FC-cuts  because  of  their  technological  significance. 
Also  given  are  charts  of  the  mode  spectra  in  the  region  of  the  thickness  modes 
and  the  modal  temperature  coefficients. 


*See  list  of  references  beginning  on  p.  39. 
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The  force  coefficient  data  are  compared  with  theoretically  predicted 
values  obtained  from  a variational  principle  applied  to  an  anisotropic  disc 
supported  at  two  diametric  points.  This  analysis  departs  from  previous 
treatments3 »5,19, 21 ,23,24  in  two  major  respects:  1)  the  isotropic  stress 

pattern  is  replaced  by  the  more  accurate  anisotropic  stress;  2)  the  elastic 
problem  is  treated  for  the  general  trlclinlc  symmetry,  rather  than  for  the 
monoclinic  symmetry  appropriate  to  rotated-Y-cuts . 

These  investigations  verify  the  predicted  superiority  of  doubly  rotated 
quartz  plates  over  the  conventional  AT-cut,  with  respect  to  the  force- 
frequency  effect,  and  provide  further  motivation  for  their  continued 
development  and  utilization. 

DOUBLY  ROTATED  CRYSTAL  PLATES 

Doubly  rotated  crystal  plates  are  the  most  general  kind  of  one- 
dimensional thickness-mode  vibrator.  The  orientation  is  uniquely  specified 
by  two  angles — <p  and  0.  Following  the  usual  convention, ^6*  the  orientation 
with  respect  to  the  crystallographic  axes  is  described  as  (YXw&)4>/6. 

Examples  of  singly  and  doubly  rotated  cuts  are  shown  in  Figure  1 along  with 
the  angles. ^ Shown  in  Figure  2 is  the  locus  of  zero  first-order  temperature 
coefficient  (ZTC)  for  quartz  resonator  plates.  In  quartz,  the  first  rotation 
lovers  the  apparent  symmetry  from  trigonal  to  monoclinic;  the  second  rotation 
further  lovers  it  to  triclinic. 

STATIC  FREQUENCY-TEMPERATURE  BEHAVIOR 

For  quartz  plates  on  the  upper  (0  > 0)  ZTC  locus,  the  static  frequency- 
temperature  (f-T)  curve  exhibits  a cubic  behavior.  The  AT-cut  is  the  classi- 
cal example.  Whereas  the  inflection  temperature  (the  temperature  half  way 
between  turnover  points)  occurs  at  room  temperature  for  the  AT-cuts,  this 
temperature  increases  steadily  with  increasing  angle  <p>  becoming  48°C  at  the 
FC-cut,  74°C  at  the  IT-cut,  95°C  at  the  SC-cut,  and  157°C  for  the  rotated-X- 
cut.23  Frequency-temperature-angle  curves  for  the  AT-  , FC-  , and  SC-cuts 
are  given  in  Figures  3,  4,  and  5,  respectively.  Typical  f-T  curves  for  the 
SC-cut  are  shown  in  Figure  6 for  the  c-  and  b-modes.  Dots  represent  measured 
points.  For  the  b-mode  (the  faster,  quasi-shear  mode,  classed  as  undesired), 
the  first-order  temperature  coefficient  is  -25.1  X10-8/K.  The  c-mode  (slower; 

quasi-shear,  desired  mode)  curve  was  fit  by  quintic  least-squares  to  yield  the 
following  coefficients  at  25°C: 

a - 0.38  X10“6/K  d - 136.  X10"15/K4 

b - -11.4  X10_9/K2  e - 175.  X10"18/k5 

c - 26.8  X10_12/K3 

The  SC-cut  is  seen  to  have  its  inflection  temperature  around  100°C,  so  that  it 
would  normally  be  operated  around  its  lower  turnover  temperature,  where  the 
upper  turnover  is  used  for  the  AT-cut.  Compared  with  the  AT-cut,  the  SC-cut 
is  also  flatter,  so  a given  temperature  control  will  correspond  to  a smaller 
frequency  deviation. 

*See  list  of  references  beginning  on  p.  39. 
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FIGURE  2.  LOCUS  OF  ZERO  TEMPERATURE  COEFFICIENT  IN  QUARTZ. 


AT-CUT 


AT  ■ T~Tj  (kelvins) 

-ANGLE  CHARACTERISTICS  OF  AT-CUT  QUARTZ  RESONATORS. 


AT“T"Tj  (kelvins) 

FIGURE  5.  F REQUEi'l CY-TEfIPERATURE -ANGLE  CHARACTERISTICS  OF  SC-CUT  QUARTZ  RESONATORS 


MODE  SPECTROGRAPHS6 

A wideband  mode  spectrograph  Is  shown  in  Figure  7 for  an  SC-cut  crystal. 
The  modes , In  order  of  frequency,  are  denoted  as  nW , where  m Is  the  mode 
type  (a,  b,  or  c)  and  M Is  the  order  of  the  harmonic.  The  sequence  shown  Is 

cjl)  b£l)  a{l)  cp)  bp)  cP)  bP)  ap)  and  c(7).  The  thlcluiess-sheai  mode 
TSi  at  cutoff  corresponds  to  the  c-mode;  the  thickness-twist  TT^  at  cutoff 
corresponds  to  mode  b , and  mode  a is  the  thickness-stretch  mode,  TE.  The 
spacings  and  amplitudes  measured  agree  closely  with  those  calculated  last 
year.6  Table  1 extends  the  calculations  to  higher  harmonics  The  ratios  of 
frequencies  are  for  unelectroded  antiresonance  frequencies. 

In  Figure  8 the  spectrum  in  the  vicinity  of  the  c-mode  Is  shown;  a 
narrowband  plot  about  the  b-mode  resonance  Is  given  in  Figure  9.  From 
Figures  7,  8,  and  9 one  sees  how  very  clean  the  spectrum  is,  even  for  the 
harmonic  modes.  The  flat  SC-plate  used  for  this  experiment  had  the  following 
measurements : 

plate  diameter  ” 14.18  mm; 

electrode  diameter  $e  “ 5.0  mm  (keyhole  pattern); 
mass  loading  (plateback)  y “ 1.8Z; 
c-mode  fundamental  frequency  ■ 5.937  MHz. 

Although  one  would  suspect  that  the  SC  orientation,  because  of  Its  lower 
synaetry,  would  have  a more  complicated  unwanted  mode  spectrum.  It  appears 
that  energy  trapping^? »^8*  can  be  applied  readily  to  these  plates,  although 
the  optimum  electrode  shaped*  and  plateback  relations  are  not  available  at 
the  present  time. 

ELECTRICAL  CHARACTERISTICS 

For  the  Ideal  case  of  an  infinite  plate  having  a uniform  distribution  of 
motion  laterally,  the  SC-cut  c-mode  yields  the  following  physical  and  electri- 
cal characteristics^  .13 . 

e Dielectric  permittivity  of  static  capacitance  CQ: 
e - CQt&/t±  - 39.8  pF/m 

a Motional  time  constant:  T ■ R C ■ P T - 11.7  fs 

l ll  ll 

a Capacitance  ratio  of  fundamental:  r'l)  ■ C0/C^  « e/I*  ■ 496 

e Frequency  constant:  N ■ 1.797  MHz-mm 

e Motional  capacitance  constant  (permittivity  of  motional  capacitance) : 

T - C t„/A  - 80.3  fF/m 

l l ■ e 


*See  list  of  references  on  p.  42. 
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TABLE  1.  SEQUENCE  AND  STRENGTHS  OF  SC-CUT  MODES 


1 > 

I ' 1 

FIGURE  8.  NARROWBAND  MODE  SPECTROGRAPH  FOR  THE  SC-CUT  c-MODE. 


b MODE 


Motional  resistance  constant  (resistivity  of  motional  resistance) : 


f ** 

R A / 1 

= 146  fl-mm. 

1 

1 e a 

In  the  above,  Ae 

= x$~/4  : 

Ls  the  electrode 

thickness. 

We  pick  for 

illustration  the  measured 

• 

» 5.937 

MHz 

. t-1 

= 3.304 

-1 

HMD 

a 

• $ 

xa 

= 14.18 

mm 

• C 

o 

= 3.79 

pF 

* C1 

= 5.91 

fF 

• R. 

= 55 

0. 

1 

From 


hese  measurements  we  deduce  the  following  effective  values: 
46.9 

(I)' 


V*. 


r' 

T 


641 

325  fs 


§e  (effective)  = (4taCo/ire) 


6.06  mm 


A / 1 
e a 


Q 


95.3  mm 

=62.0  fF/ra 
= 5.24  fl-m 

= 1/u  r = 82,500 

= T (effective) /r  (theoretical) 
1 1 


77.3%. 


In  the  above  the  effective  electrode  diameter  ♦ was  determined  from  the 

0 

measured  value  of  C . The  measurement  of  Cq  and  was  carried  out  on  an  RX 
Meter  (schering  bridge)  using  a frequency  synthesizer.  In  the  vicinity  of 
resonance,  the  parallel  capacitance  versus  frequency  curve  approximates  a 
hyperbola  very  closely,  and  C and  are  obtained  at  the  same  time.  Our 
measurements  were  fit  by  leas?-squares  to  improve  the  fit.  The  presence  of 
stray  capacitance  effects  was  also  carefully  minimized.  The  nominal  elec- 
trode diameter  was  5 mm,  whereas  the  effective  value  is  about  6 ram;  the 
discrepancy  may  be  partly  due  to  the  electrode  tabs,  with  some  component  due 
to  stray. 


14 


It  is  interesting  to  compute  the  values  of  Cp  R^,  and  Q that  would 
result  in  the  ideal  case  of  uniform  distribution  of  motion  and  all  losses  due 
solely  to  bulk  wave  attenuation  due  to  the  viscosity.  The  results  are: 

• = CQ/r^  (theoretical)  - 7.64  fF 

• R,  m P , (theoretical)  • t /k  * 1.53  ft 

11  a e 

• Q ■ l/o^  t (theoretical)  = 2.29  XIO^. 

Taking  into  account  the  nonuniform  distribution  of  motion  with  T,  the 
motional  resistance  would  become  (again  assuming  all  loss  is  due  to 
viscosity) : 

• (theoretical)  /'V  = 1.98  Cl. 

Most  of  the  loss  in  this  resonator  is  due  to  factors  other  than  the  bulk  wave 
viscosity.  In  fact,  since 

• t (ef fective) /t  (theoretical)  ■ 27.8, 

1 1 

the  bulk  wave  component  is  only  (1/26.8)  or  less  than  4%  of  the  total. 

At  present,  the  apportionment  of  the  losses  in  doubly  rotated  resonators 
is  largely  an  open  question.  It  was  observed,  however,  that  those  SC-cut 
resonators  having  the  electrode  tabs  running  along  Z"  had  resistance  values 
from  three  to  ten  times  smaller  than  those  with  tabs  aligned  along  X".  This 
points  to  energy  being  carried  away  to  the  mounting  at  the  plate  edges. 
Associated  with  this  is  the  fact  that  the  doubly  rotated  cuts  in  quartz  do 
not  support  pure  modes,  such  as  the  AT-cut  mode.  This  means  that  the  SC-cut 
c-mode  consists  in  particle  motion  that  is  out  of  the  plane  of  the  plate. 
Considering  the  plate  axes  as  the  reference  set,  and  performing  4 and  0 
rotations  with  respect  to  this  plate  reference  coordinate  system,  the  angles 
and  of  the  particle  displacement  for  the  SC-cut  c-mode  are: 

• <J> . - 8.57° 

a 

• e * 13.7° 

a 

The  out-of-plane  displacements  are  not  negligible  and  certainly  represent, 
for  plano-plano  resonators  as  described  here,  a significant  mechanism  for 
energy  loss  to  the  ambient  fluid  and  the  mounting  system.  Beveling  or  con- 
touring the  resonators  and  operation  in  a vacuum  enclosure  would  decrease  the 
size  of  these  losses. 

Energy  trapping  is  another  useful  tool  for  keeping  the  Q high,  as  well 
as  maintaining  a clean  mode  spectrum;  the  energy-trapping  rules  for  doubly 
rotated  plates  have  yet  to  be  worked  out. 

Additional  work  must  also  be  done  to  calculate  the  frequency-wavelength 
dispersion  diagrams  for  doubly  rotated  resonators.  If  the  ratio  ♦ It  used 
happens  to  fall  at  a position  where  there  is,  e.g.,  strong  coupling  to  flexure 
waves  propagating  laterally  across  the  plate,  then  energy  would  be  carried 
from  the  desired  mode  into  the  unwanted  motion. 
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MATHEMATICAL  MODELING 


An  overview  of  the  present  theoretical  work,  compared  to  past  treatments, 
is  shown  in  Figure  10. 

Past  theoretical  analyses3,5,19,21  of  static  mechanical  stress  bias 
effects,  in  general,  and  of  force  sensitivities,  in  particular,  have  been  two- 
step  calculations.  First,  linear  elastic  solutions  for  the  distribution  of 
static  mechanical  stress  bias  in  the  resonator  blank  were  obtained  assuming 
that  a-quartz  is  isotropic.  These  isotropic  static  solutions  for  the  stress 
at  the  blank  center  were  then  used  to  calculate  resonant  frequency  shifts  in 
nonlinear  wave  propagation  calculations  which  included  the  correct  anisotropy 
of  quartz  and  third-order  elastic  constant  effects. 

The  present  work  is  a two-step  calculation  where  an  attempt  is  made  for 
a better  solution  to  the  static  problem.  The  calculus  of  variations  (essen- 
tially the  Rayleigh-Ritz  method) 50*  is  U8ed  to  find  approximate  solutions  to 
the  anisotropic  static  stress  problem.  The  approximate  static  solution  for 
the  stress  at  the  blank  center  is  used  in  a nonlinear  elastic  wave  propagation 
code  to  calculate  resonant  frequency  shifts. 

The  present  theoretical  results  for  the  singly-rotated  AT-  and  BT-cuts 
provide  a much  improved  comparison  with  published  experimental  results  than 
the  earlier  theoretical  results.  The  theory  is  used  to  calculate  force  sensi- 
tivity coefficients  for  the  doubly  rotated  family  containing  the  important 
AT-  , FC-  , IT-  , and  SC-cuts.  The  results  provide  the  crystal  designer  with 
the  appropriate  azimuthal  angle  to  mount  the  resonator  on  a two-point  mount 
for  minimum  force  sensitivity. 


One  calculational  approach  that  has  been  useful  on  numerous  occasions  for 
calculating  stress  patterns  in  static  and  vibrating  elastic  material  bodies  is 
the  calculus  of  variations .51,52*  The  method  is  approximate,  although  the 
closeness  to  which  the  approximate  solution  can  be  brought  to  the  actual  solu- 
tion is  a matter  of  degree  depending  on  the  choice  of  trial  (basis)  functions, 
available  computer  size,  and  patience.  The  method  for  static  problems  amounts 
to  formulating  the  total  elastic  stress-strain  energy  stored  in  a given  body 
for  the  given  boundary  conditions  and  trial  functions  and  then  adjusting  the 
trial  functions  to  minimize  the  stored  elastic  energy.  The  approximations 
obtained  with  the  method  are  somewhat  better  for  the  elastic  energy  values 
than  for  the  stress  distributions,  but  sufficient  accuracy  of  the  stress  dis- 
tributions can  be  obtained  for  practical  considerations. 


The  total  stored  elastic  energy  L is  given  by 


L ■ i jffifav  - ff*2  Fiui 

V s 

Here  V and  S are  the  resonator  blank  volume  and  surface,  Cy^  is  the  elastic 
stiffness  tensor  in  engineering  notation,  Fj  is  the  distribution  of  force  per 
unit  area  acting  on  the  surface,  Uj  is  the  elastic  displacement  vector.  We 
use  a cartesian  coordinate  system  for  the  plate,  and 

*See  list  of  references  on  p.  42. 
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MATHEMATICAL  MODELING  OF  FORCE-FREQUENCY  EFFECT 


Prior 


Present 


FIGURE  10. 


• Isotropic  initial  stress  field  assumed. 

• Analysis  limited  to  rotated-Y-cuts  (YX£)9. 

• Anisotropic  initial  stress  field  obtained 

BY  CALCULUS  OF  VARIATIONS. 

• Solution  applied  to  doubly  rotated  cuts 

(YXux0$/0. 


SYNOPSIS  OF  FORCE-FREQUENCY  EFFECT  MODELING. 
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(2) 


’ij 


<U 


i.j 


Uj,i>/2‘ 


Here  X and  y run  1-6,  and  i and  j run  1-3.  is  related  to  by  the 
conventional^!  relations  between  engineering  and  tensor  notations. 


The  present  calculations  treat  a circular  resonator  blank  of  diameter  d 
and  thickness  x described  by  the  IRE  standard^**  notation  (YXw£)$/0.  We  con- 
sider only  contoured  or  energy-trapped  resonator  designs  where  the  vibra- 
tional acoustic  energy  is  restricted  to  the  vicinity  of  the  blank  center. 
Contouring  effects  are  ignored  for  the  static  stress  distribution  calculation, 
however,  so  that  the  much  simpler  problem  of  a flat  circular  plate  can  be 
solved  for  the  static  stress  distribution.  This  simplification  is  warranted 
because  the  thin  resonator  blanks  used  for  thickness  shear  resonators  allow 
the  assumption  that  the  thickness  dimension  is  small  enough  for  a plane 
stress  problem. 


In  the  case  of  plane  stress,  X and  y run  1,  3,  5 and  i and  j run  1,  3, 
in  Equations  1 and  2 (X2  is  blank  thickness  direction,  x^  is  l,  X3  is  w in 
the  standard  notation).  Also,  is  replaced  by  y.  , the  planar  elastic 
stiffness  coefficients  expressed  in  the  plate  coordinate  system.  Hence: 


■ 3®  Wui  • JJ “ Fiui 


(3) 


The  variational  method  involves  substituting  a linear  superposition  of 
trial  functions  for  into  Equations  2 and  3,  carrying  out  the  integrals  in 
Equation  3,  and  minimizing  the  resulting  expression  (differentiating  the 
expression  with  respect  to  a given  coefficient  and  setting  that  equal  to 
zero)  with  respect  to  the  coefficients  of  the  trial  functions.  The  problem 
then  becomes  a linear  algebra  problem  in  the  coefficients.  The  choice  of 
trial  functions  must  be  such  that  they  represent  a pointwise  complete  set 
over  V and  S.  If  the  trial  functions  already  satisfy  some  aspect  of  the 
problem  such  as  the  differential  equation  or  boundary  conditions,  the  number 
of  trial  functions  needed  for  adequate  convergence  is  small.  For  the  present 
case,  we  take  the  easy  way  out  and  choose  a simple  power  series  expansion  for 
and  rely  on  the  computer  to  handle  large  numbers  of  trial  functions.  The 
trial  functions  chosen  for  the  two-point  problem  are 
MNPQ 


Ul  = 


E 


. v2n+l  2n 
A x x + b 
mn  1 3 


2p  ^qfl 

pq  1 3 


(4) 


m,n,p,q=0 


and 


U„ 


x^1  x 2q  + n x. 
1 3 tu  1 


2t 


*rj 


(5) 


The  value  for  is  set  equal  to 

Ft  - F/(6t), 
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(6) 


where  F is  the  Inwardly  acting  force  applied  to  the  opposite  ends  of  a blank 
diameter,  t * ta  is  the  plate  thickness,  and  6 is  some  length  dimension  which 
is  small  relative  to  the  blank  circumference  (point  force).  All  the  inte- 
grals can  be  carried  out  in  Equation  3,  the  surface  integral  being  treated  in 
the  limit  of  a point  force  (6  -*•  0) . After  differentiation  with  respect  to 
A^,  B , and  Dtu,  the  resulting  linear  algebra  problem  has  a solution  for 
Ajjjj,  , and  Dtu,  which  scales  with  F.  Thus,  the  static  stress  distribution 

can  bePsolved  for  any  known  y A computer  code  for  the  linear  algebra 
problem  was  written  which  included  rotation  of  the  quartz  elastic  tensor  to 
obtain  y.  and  arbitrary  selection  of  M,  N,  P,  Q,  T,  and  U.  From  symmetry 
arguments;  one  is  led  to  use  groups  of  the  trial  functions,  adding  new  groups 
until  satisfactory  convergence  is  obtained.  The  groups  are  defined  by  (0,R) , 
where  the  group  includes  all  pairs  of  (m,n) , (p,q),  and  (t,u)  with  the  first 
member  increasing  from  zero  to  R in  steps  of  one  while  the  second  member 
decreases  from  R towards  zero  in  steps  of  one;  e.g.,  (0,3),  (1,2),  (2,1), 
(3,0)  make  up  the  family  (0,3). 


The  resulting  solution  for  the  static^stress  at  the  plate  center  using 
the  published  constants  for  quartz^  f8  incorporated  into  a previously 

described  computer  code5»5  which  calculates  the  resonant  frequency  shift 
caused  by  elastic  nonlinearities  (third-order  elastic  constant  effects). 54 
We  use  here  the  definition  for  the  force  sensitivity  coefficient  Kf  defined 
as^O 

Af  FN0 

T " Kf  dT  <7> 

where  N0  is  the  frequency  constant.  The  plate  diameter  is  d ■ *a.  Units  are 
Kf  in  m*sec/N,  N0  in  m/sec,  d and  x in  m,  and  f in  sec'*  or  Hz.  Kf  is  posi- 
tive if  frequency  increases  upon  application  of  a compressive  force. 


The  direction  of  the  applied  forces  F is  important  because  of  the  aniso- 
tropy of  the  quartz  nonlinear  elastic  problem  and  because  of  the  quartz 
anisotropy  in  the  static  stress  problem  (the  latter  is  ignored  in  earlier 
isotropic  static  solutions) . We  choose  to  follow  the  earlier  experimental 
work  by  presenting  results  for  the  Kf  of  a given  (YXwi)$/0  cut  as  a function 
of  azimuthal  angle  ty.  The  azimuthal  angle  is  measured  in  right-hand  conven- 
tion about  x'2  from  xj  (or  A)  in  the  plane  of  the  plate:  positive  angle  going 

from  xj[  to  -x-j*  as  seen  in  Figure  11.  The  direction  of  the  applied  force  is 
therefore  specified  by  Xj£  when  it  undergoes  the  rotations  (YXwlt)$/0/4>.  The 
angles  $,  0,  and  tji  are  just  the  euler  angles  encountered  in  classical 

mechanics. 


B.  Results 

The  program  was  first  tested  for  isotropic  blanks.  The  solution  at  the 
center  of  the  plate  converges  rapidly  (four  significant  figures)  to  that 
found  analytically, 56  i.e.,  6F/irdx  compression  along  the  diameter  aligned 
with  (inwardly  directed)  F and  2F/irdx  tension  along  the  diameter  perpen- 
dicular to  F. 

Calculations  were  carried  out  for  the  AT-  and  BT-cuts  since  there  exist 
large  amounts  of  force  sensitivity  data  and  considerable  discrepancies 
between  the  data  and  earlier  theoretical  results  for  these  cuts.  The  results 


*See  list  of  references  on  p.  43. 
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for  the  zero  temperature  coefficient  thickness  shear  c-mode  are  shown  in 
Figure  12  as  Kf  vs.  tfi.  The  "isotropic  assumption"  in  Figure  12  is  the  result 
obtained  in  the  present  computer  codes  when  the  isotropic  solution  (-6F/irdx, 
2F/irdr)  is  used  for  the  static  stress  pattern  at  the  plate  center.  This 
result,  which  is  representative  of  the  best  results  obtained  before  the 
present  work,  is  numerically  equivalent  to  the  results  of  Lee,  et  al.^l  Note 
the  discrepancies  between  the  isotropic  assumption  results  and  the  experimen- 
tal summary  provided  by  Ratajski;20  namely,  the  difference  in  value  for  Kf 
at  <l>  = 0 for  the  AT-cut,  and  the  complete  failure  to  predict  the  dip  in  the 
BT-cut  results.  As  seen  in  Figure  12,  the  present  calculations  using  the 
variational  treatment  are  quantitatively  accurate  at  iji  = 0 for  the  AT-cut  and 
predict  a dip  for  the  BT-cut.  The  calculations  leading  to  Figure  12  involved 
using  all  the  families  up  to  and  including  (0,5),  making  21  values  each  for 
Amn»  Bpq , and  Dtu  for  a 63  x 63  linear  algebra  problem.  The  addition  of  the 
(0,5) family  only  changed  the  numerical  answer  for  the  stress  at  the  plate 
center  in  the  third  significant  figure,  so  that  some  idea  of  convergence 
would  be  provided. 

Figure  13  contains  present  results  for  the  thickness-shear  (mode  c) 

IT-cut  and  previously  published  experimental  results. ^ The  qualitative  fea- 
tures of  the  experimental  data  are  reproduced  by  both  the  isotropic  assumption 
results  and  the  full  variational  result.  In  view  of  the  results  in  Figure  12, 
where  more  experimental  data  exist,  we  consider  the  variational  result  to  be 
the  more  accurate.  The  quantitative  discrepancy  between  experiment  and  theory 
is  not  as  serious  as  might  appear  in  Figure  13  because  there  is  a factor-of- 
two  scale  expansion  between  Figures  11  and  13  to  account  for  the  generally 
smaller  Kf  values  of  doubly  rotated  cuts. 

A series  of  calculations  has  been  carried  out  for  the  zero  temperature 
coefficient  branch  (mode  c)  subset  within  (YXwi)$/0:  4>  = 0°  (AT-cut); 

<J>  =*  10°;  $ = 15°  (FC-cut);  <f>  = 19.1°  (IT-cut);  <(>  * 21.9°  (SC-cut);  * - 30° 
(rotated-X-cut) . The  results  for  Kf  vs.  \|i  are  shown  in  Figure  14.  The  curves 
drawn  in  Figure  14  are  presented  separately  in  Figures  15-20. 

The  family  of  curves  for  Kf  vs.  ip  in  Figure  14  are  the  general  theoreti- 
cal results  for  this  report.  The  curves  should  be  useful  in  choosing  the 
optimum  location  of  bonding  pads  for  two-point  mounts  to  minimize  the  effects 
of  mounting  stresses  and  some  effects  of  shock  and  vibration. 

In  Table  2 are  given  the  computed  angles  at  which  Kf  is  zero,  as  func- 
tion of  $.  Table  3 lists  the  computed  locations  and  values  of  the  extrema  of 
Kf  as  function  of  <p . 


TABLE  2.  ZEROS  OF  Kf 


r 

<J>° 

r 

0 

64.7 

, 115.3 

19.1 

79.3 

, 163.1‘ 

10 

68.5 

, 125.2 

21.9 

81.6 

, 171.9 

15 

74.8 

, 148.8 

30 

79.3 

, 184.3 

21 


ISOTROPIC  ASSUMPTION 

VARIATIONAL  RESULT 

RATAJSKI  SUMMARY  (EXPERIMENTAL) 


AZIMUTHAL  ANGLE  (DEGREES) 

FIGURE  12.  FORCE-FREQUENCY  RESULTS  FOR  AT-  AND  BT-CUT  QUARTZ. 


FREQUENCY  SHIR  COEFFICIENT  (10'A:>  M.SEC/N) 


FIGURE  13,  FORCE-FREQUENCY  RESULTS  FOR  IT-CUT  QUARTZ. 


<LO  m-a/N 


*<L o »v) . s/N 


TABLE  3.  LOCATION  AND  VALUES  OF  Kf  EXTREMA 


*° 

r 

Kf  (max) 

r 

Kf  (min) 

0 

0 

24.5 

90 

-11.5 

10 

17.6 

21.4 

93.2 

-10.8 

15 

26.0 

18.2 

98.8 

- 8.9 

19.1 

36.6 

15.3 

118.1 

-10.1 

21.9 

44.3 

14.7 

131.6 

-14.3 

30 

55.6 

14.7 

130.0 

-23.9 

(Kf  in  10-15  m*s/N) 

EXPERIMENTAL  CONFIRMATION 

Apart  from  the  IT-cut  results  obtained  in  I960,*  no  experimental  results 
were  available  for  comparison  with  the  theoretically  predicted  force- 
frequency  curves  of  doubly  rotated  cuts  of  quartz,  shown  in  Figure  14.  In 
order  to  obtain  these  results,  the  experimental  apparatus  of  Figure  21  was 
constructed.  It  contains  provisions  for  the  precise  rotation  of  the  angle  iji 
by  means  of  a vacuum  chuck  for  holding  the  crystal  (and  maintaining  strict 
crystal  vertical  alignment)  and  high-ratio  reduction  gear.  Five  mil  wires 
are  bonded  to  the  crystal  and  serve  for  the  electrical  connections.  To  the 
left  in  the  figure,  the  connections  are  brought  out  to  a Crystal  Impedance 
Meter  (RFL  model  459  ” TS-33owith  low  drive  modification^?*) . The  frequency 
is  recorded  from  a counter. 

Micrometer  adjustments  are  available  for  assuring  accurate  alignment  of 
the  various  portions  of  the  jig.  Force  application  was  made  by  a movable  rod 
on  which  calibrated  masses  were  applied;  the  mass  of  the  rod  was  taken  into 
account.  The  overall  setup  is  shown  in  Figure  22. 

The  crystals  considered  here  had  orientations  (YXwZ)$/0,  with  $ ■ 10#, 
15°,  and  21.9°,  and  (p  such  that  the  units  had  zero  temperature  coefficients 
(<t>  = +34°) . All  crystal  units  were  provisionally  scored  near  the  X?  axis 
(projection  on  the  Z,  or  optic,  axis).  The  true  location  of  X'j'  with  respect 
to  the  score  mark  was  later  determined  for  each  crystal  to  within  1°  by  a 
conoscope,  and  enabled  to  be  accurately  known.  (This  procedure  was  also 
applied  to  the  still-extant  IT-cut  crystal  that  was  described  in  Reference  1. 
The  new  measurement  disclosed  that  the  true  X'j  axis  as  seen  in  the  conoscope 
was  a full  19°  in  error  with  respect  to  the  score  mark  on  the  crystal.  Thus, 
the  curve  in  Figure  18  of  Reference  1 ought  to  be  translated  to  the  right 
with  respect  to  the  graph  axes,  so  that  the  zeros  occur  at  85°  and  163°. 

This  finding  was  a welcome  resolution  of  a disturbing  discrepancy  between 
theory  and  experiment!) 

Measurements  were  made  in  \p  intervals  of  10°.  Three  readings  were  taken 
at  each  \|i.  First,  the  frequency  of  the  unloaded  vibrator  was  measured.  Then 
the  weight  was  lowered  gently,  and  the  loaded  frequency  was  recorded. 

Finally,  the  frequency  with  the  weight  removed  was  measured.  The  first  and 
third  frequencies  were,  in  all  cases,  at  most  one  or  two  Hz  apart.  Kf  was 

”&ee  list  of  references  on  p.  43. 
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then  calculated  from  equation  (7)  using  th®  thickness , force**  applied, 

change  Af  with  applied  weight,  Plate  ^^®  cut8  considered,  the  anti- 

and  crystal  frequency  constant  N0.  For  th- 
resonance  frequency  constants  are 


10°, 

15“, 

21.9°, 


1690  m/s 

1726 

1797 


and  25  show  the  experimental  curves  for  the  c-mode  cuts 
Figures  23,  24  and  25  Included  in  each  figure  is  a curve, 

the  theoretical  result  obtained  from  the  variational  proce- 

trJ’JZrt  represent  data  •*«».  encountered. 

Date  acattar  la  worst  for  the  » - 10"  units, 
explanation  for  this  cannot  at  present  vary  wlth  ^ indicating, 

observed  that  the  unloaded  ftecuenc  e pp  proximities  of  the  con- 

e.g.,  . stray  capacitance  effect  do.  Mttaoh  angl^g  „„„geTO»t  to 

necting  wires  to  the  fixture.  curve  shown  in  Figure  25  without  the 

reduce  this  effect  produced  e f modification  improves  the 

error  bars.  As  may  be  seen  in  the  ^he  th^Stical  In  Figures  23,  24, 

agreement  of  the  experiment  curve  theory  is  generally  quite  good, 

and  25,  the  agreement  between  VSnd.  The  overall 

especially  when  the  smallness  do  the  general  features  in 

magnitudes  predicted  and  ^ observed  in  the  AT-  and  BT-cut 

each  case.  The  symmetry  about  \p  aWnt  in  the  doubly  rotated  cut  curves 
curves  (of.  Figure  12)  is  soon  to  be  absent  in  the  doub  y of  try, 

SK  in ' the  * plane^f  “^uT^s,  1.  oSt  of  the  plane  for  cuts 

having  <fr  I*  0°. 

_j  oc  gjtth  Fiaure  12  reveals  that  the  peak- 

to-peakCexcursions°of  with  * for  the  SesHutfto 

“zs-jzr-tsz'zzsx  " “ « - » “ 

that  the  average  force  coefficient 

l rn  (8) 

<Kf>“  KfW 

for  the  FC-cut 

as  the  ZTC  orientation  for  which  planar 

stress  produces  no  frequency  change. 

I..  a nroUmlnarv  set  of  measurements  of  the 

Figure  26  gives  the  r®sul^s  ® SC-cut  7 This  mode  has  antiresonance 

force  sensitivity  for  the  b-mode  of  the  SC-cut. 

frequency  constant 

Nq  * 1977  m/s, 

and  a temperature  co.fflcl.nt  of  approximately  -25  or  -26  ppm/g.  Because  of 
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<H0* 


AZIMUTHAL  AKCLEf DECREES) 

FIGURE  23.  EXPERIMENTAL  AND  THEORETICAL  RESULTS  FOR  phi  EQUALS  10  DEGREES. 


(FC-cut) 


AZIMUTHAL  ANCLE (DECREES) 

FIGURE  2A,  EXPERIMENTAL  AND  THEORETICAL  RESULTS  FOR  phi  EQUALS  15  DEGREES. 


AZIMUTHAL  ANGLE  ( DECREES) 

FIGURE  25,  EXPERIMENTAL  AND  THEORETICAL  RESULTS  FOR  phi  EQUALS  21.9  DEGREES. 


the  large  temperature  coefficient,  the  force  coefficient  measurements  are  not 
as  easy  to  make  as  for  the  c-mode.  The  results  reported  In  Figure  26  should 
be  regarded  as  tentative;  they  disclose,  however,  that  Kf  exhibits  no  zeros 
along  the  ip  axis. 

CONCLUSIONS 

Doubly  rotated  quartz  resonators  have  been  considered  as  to  frequency- 
temperature  behavior,  mode  spectrum  content,  and  primarily  with  respect  to 
the  force-frequency  effect.  It  was  found  that  the  f-T  curve  shifts  upward 
In  temperature  with  increasing  <P  angle  and  becomes  flatter.  The  mode  spec- 
trum is  complicated  by  the  presence  of  all  three  thickness  modes,  but  their 
enharmonic  overtones  may  be  very  adequately  suppressed  by  means  of  energy 
trapping. 

In  the  force-frequency  portion,  the  calculus  of  variations  has  been 
used  to  obtain  solutions  for  the  static  stress  distributions  caused  in 
circular  quartz  resonator  blanks  by  forces  acting  along  a diameter.  The 
resulting  static  stress  distributions  have  been  used  to  predict  the  force 
sensitivity  of  the  resonant  frequency  of  thickness  shear  resonators.  The 
present  results  are  shown  to  agree  more  favorably  with  experiment  than 
results  obtained  by  previous  theories  where  the  quartz  was  assumed  to  be 
isotropic  for  the  static  stress  solution.  General  results  are  given  by  the 
doubly- rotated  family  (YXv£)4>/0  along  the  technologically  Important  zero 
temperature  branch  which  contains  the  AT-  , FC-  , IT-  , and  SC-cuts.  The 
results  will  be  useful  for  choosing  the  location  of  bonding  pads  for  two- 
point  mounts  so  that  mounting  stress  effects  can  be  minimized. 
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APPENDIX  A 


RESONATOR  PREPARATION 

The  resonator  plates  used  for  these  Investigations  were  of  three 
orientations : 

4>  - 10°  , 0 * 34°50* 

15°  34°16'  (FC-cut) 

21.9°  33°54'  (SC-cut). 

Each  orientation  group  consisted  of  four  crystals.  Initial  thicknesses 
varied  from  ta  - 0.635  nan  to  1.397  mm.  All  crystal  blanks  were  rough-lapped 
to  a thickness  of  0.559  mm  with  #400  mesh  silicon  carbide  in  a water  slurry. 
The  plates  were  then  reduced  to  the  operating  thickness  according  to  the 
following  schedule: 

Grit  Size 


# 600 
12  pm 
5 pm 
3 pm 

Polishing 
Compound 
(cerium 
oxide) 

The  final  plate  thickness  brings  the  c-mode  fundamental  resonance  into  the 
neighborhood  of  6 MHz. 

In  the  lapping  process,  the  crystals  were  grouped  according  to  orienta- 
tion and  thicknesses  to  within  +0.013  mm  for  each  of  the  above  stages.  These 
groups.  In  turn,  were  cemented  witi.  a wax-rosin  compound  to  a 75  mm  diameter 
optical  flat  (flatness  better  than  /10)  during  lapping  as  a means  of 
assuring  orientational  integrity  and  parallelism. 

All  blanks  were  rounded  after  lapping.  The  plates  were  loafed  together 
in  standard  fashion,  and  the  diameter  reduced  to  14.18  mm  for  the  SC-cut 
units  and  13.97  mm  for  the  <p  ■ 10°  and  15°  cuts.  The  final  diameter  was 
reached  when,  upon  microscopic  inspection,  all  hairline  cracks,  chips,  and 
scratches  were  removed. 

The  cleaning  procedure  consisted  in  washing  the  blanks  in  detergent  and 
rinsing  them  under  running  tap  water,  as  a first  step.  (Cleaning  was  limited 
to  this  first  step  in  the  washings  that  took  place  between  successive  stages 
of  the  lapping  procedure  discussed  above,  and  during  each  reversal  of  the 
plates  on  the  optical  flat.) 

Cleaning  of  the  polished  resonator  blanks  prior  to  electrode  plating 
consisted  of  the  further  steps: 


Material  Removed 
0.076  mm 
0.076  mm 
0.051  mm 
0.038  mm 
0.013  mm 


Plate  Thickness 

0.483  mm 
0.407  mm 
0.356  mm 
0.318  mm 
0.305  mm 


Ho** 
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• Inserting  blanks  into  specially  constructed  stainless  steel  wire  jigs. 

• Degreasing  in  boiling  acetone  or  methanol  for  3 minutes. 

• Rinsing  In  hot  distilled  water  in  two  successive  stages. 

• Agitating  in  ultrasonic  bath  with  hot  ethyl  alcohol  for  1 or 
2 minutes. 

• Blowing  dry  with  hot  air  (blanks  still  in  jig). 


After  the  cleaning,  the  blanks  were  removed  from  the  holder  jig  by 
teflon-coated  tweezer  and  placed  in  a UV/ozone  cleaning  facility.  This  has 
been  shown  to  be  an  effective  way  of  removing  a variety  of  contaminants  from 
quartz  surfaces.**^ 


The  blanks  were  electroded  with  argon-sputtered  gold  films  to  a diameter 
of  5 mm  with  a connecting  tab  in  a keyhole  pattern.  Plateback  values  were  as 
f ollows : 

\i  * 0.82  , <J>  - 10° 

1.32  15° 

1.82  21.9° 


Excitation  of  these  plates  was  made  possible  by  bonding  0.127  om  wires 
directly  to  the  electrode  tabs  with  electrically-conducting  cement. 
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